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We study low-temperature transport through carbon nanotube quantum dots in the Coulomb 
blockade regime coupled to niobium-based superconducting leads. We observe pronounced con- 
ductance peaks at finite source-drain bias, which we ascribe to elastic and inelastic cotunneling 
processes enhanced by the coherence peaks in the density of states (DOS) of the superconducting 
leads. The inelastic cotunneling thresholds display a marked dependence on gate voltage caused by 
different tunneling-renormalizations of the two subbands in the nanotube. Finally, we discuss the 
gate-dependent sub-gap structure observed in a strongly coupled device with odd electron occupa- 
tion. 

PACS numbers: 73.21.La, 73.23.Hk, 73.63.Fg, 74.50.+r 



Superconducting electrodes provide a useful means of 
sharpening the spectroscopic features observed in tunnel- 
ing experiments. In the superconducting phase, an other- 
wise nearly constant DOS acquires a gap of width 2A cen- 
tered at the Fermi level and characteristic sharp coher- 
ence peaks at the gap-edges =bA. These peaks transform 
a featureless metallic electrode into a high-resolution tun- 
neling probe. This widely used investigative tool^ was 
demonstrated already by Giaever's seminal work^ from 
1960 and more recently used to obtain a high resolution 
bias-spectrum of the levels in a metallic Al-nanoparticle^. 

Here we report low-temperature transport measure- 
ments in which this type of BCS-focusing promotes an 
otherwise featuresless elastic cotunneling conductance 
to sharp peaks at bias voltages Vsd = ±2A/e, corre- 
sponding to the onset of quasiparticle cotunneling. In 
the same way, inelastic cotunneling processes involv- 
ing transitions between two subbands in the nanotube 
are revealed as sharp peaks rather than steps or cusps 
in the nonlinear conductance. This sharpening of co- 
tunneling lines inside the Coulomb diamonds allows us 
to investigate more closely the tunneling-induced gate 
voltage dependence of the orbital splittinj^. Finally, 
we discuss an unusual sub-gap structure observed in 
a particularly well coupled device, signalling the im- 
portance of both multiple Andreev reflections (MAR) 
and dynamically generated bound states in spinful dots. 

A number of experiments have already investigated in- 
teresting aspects of quantum dots with superconduct- 
ing electrodes, such as supercurrenlP^, MAlP, and 
effec ts of size and charge quantiz ation in the Fabry- 
Perotl 9 | iQ | ii | i 2 | Kondd ^^ ' ^ ^ ' i5 ' i^ E3^ and Coulomb block- 
ade regimeJ^^miMiniSl^ present measurements 
performed on two different single walled carbon nan- 
otube (SWCNT) quantum dots coupled to niobium- 
based superconducting leads^A Our first device (de- 
vice A) is poorly coupled to the leads and shows reg- 
ular Coulomb blockade diamonds with clear onset of 
quasiparticle (elastic cotunneling) current. The sec- 
ond device (device B) is more strongly coupled to the 
leads and displays a four-fold degenerate shell struc- 
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FIG. 1: (Color online) Device A. (a) Bias spectroscopy at 
T = 0.3K of a SWCNT coupled to Ti/Nb/Ti leads, (b) 
Conductance within the dashed rectangle in (a). The onset 
of quasiparticle cotunneling is clearly seen at eVsd = =b2A 
(green arrows), (c) Schematic energy diagram for a viable 
cotunneling process in the center of a diamond, (d) Bias cuts 
at positions (I, II) in (b) showing conductance peaks due to 
quasiparticle (QP) cotunneling, and sequential tunneling into 
ground (GS) and excited states (ES), respectively. 



ture with both elastic and inelastic cotunneling lines 
which are sharpened by the superconducting leads. 

The devices are made by growing high quality 
SWCNT by chemical vapor deposition from predefined 
catalyst island^^^ES]^ Contacts are made of Nb-based tri- 
layers, x/Nb/x (about 5/60/10 nm), with x=Pd,Ti where 
X is thermally evaporated and the Nb is deposited by 
sputtering technique. The superconducting trilayers are 
tested via four terminal devices on the same chip and 
show transition temperatures close to T = 9K. Never- 
theless, the actual gap at the nanotube indicates a crit- 
ical temperature of Tc ~ 1.7 K, which might be related 
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to the formation of NbO (T^ - 1.4 K)P or contamina- 
tion of the lower Nb/x interface. Samples from three 
different processing rounds revealed similarly reduced Tc 
in transport measurements. However, the high qual- 
ity of the presented measurements is a promising first 
step towards Nb based SWCNT Josephson junctions. 

Figure [ija) shows the conductance versus gate, and 
bias voltage (bias spectroscopy plot) at T = 0.3 K for 
device A consisting of a carbon nanotube quantum dot 
coupled weakly to Ti/Nb/Ti leads. It reveals more than 
80 regular Coulomb diamonds, illustrating that only one 
quantum dot is defined in this high quality carbon nan- 
otube. The charging energy, [7^5 — 6meV is esti- 
mated from the height of the diamonds and no clear shell 
structure is observed^^. A region of highly suppressed 
conductance around zero bias is clearly observed for all 
gate voltages reflecting the superconducting energy gap 
of the leads. Figure [TJb) shows the Coulomb diamonds in 
the dashed rectangle of Fig. [ija). In Coulomb blockade 
(inside the diamonds) the onset of quasiparticle tunnel- 
ing is seen as horizontal lines (conductance ridges) at 
^ysd = i2A ^ ±0.55 meV (green arrows). Higher order 
Andreev reflections, which would lead to current below 
the gap, are strongly suppressed due to the poor coupling 
to the leads. Inside the diamonds, the onset of quasiparti- 
cle cotunneling corresponds to an alignment of the super- 
conducting DOS peaks and involves elastic cotunneling 
processes as depicted in Fig.jljc). At the charge degener- 
acy points, conductance inside the gap is due to Andreev 
reflections^^. Figure [ijd) shows a bias cut slightly off 
resonance (I) and further off resonance (II). The onset of 
quasiparticle tunneling at (I) involves a sequential tun- 
neling process and is therefore much stronger than in (II) . 
Peaks at higher bias are due to sequential tunneling to 
ground and excited states, respectively. 

Figure |2|a) shows a bias spectroscopy plot for the sec- 
ond, better coupled, device (device B) at T = 6.5 K 
well above the transition temperature (Tc ~ 1.7 K) of 
the superconducting Pd/Nb/Pd layer. A regular pat- 
tern is seen with three consecutive small Coulomb di- 
amonds followed by a larger diamond, reflecting the 
filling of shells consisting of two nearly degenerate or- 
bitals. Numbers in Fig. [2] denote the additional elec- 
tron number on the SWCNT for filled shells counted 
from Vgate ~ — lOV. A charging energy oi U ^ 
12meV and a level spacing of ^ 6meV are found 
from the plot. Figure [2|b) shows the conductance 
at T = 0.3 K, i.e. below T^. The lines at Vsd ^ 
zb0.55mV (green horizontal arrows) are caused by elas- 
tic quasiparticle cotunneling as illustrated by Fig. [lie). 

Inelastic cotunneling lines [Fig. |2|b) , blue arrows] are 
observed at higher Vsd foi" electron numbers N-\-l, N-\-2^ 

+ 3, but not in the fuh sheh (TV + 4 electrons). These 
lines have a marked gate voltage dependence which re- 
sembles the 'double-headed arrow' structure pointed out 
in Ref. [4l With the enhanced spectroscopy offered by 
the superconducting leads we also observe an additional 
third 'arrowhead' outside the strong elastic cotunneling 
lines appearing in every fourth (N+3) diamond, i.e. a 




FIG. 2: (Color online) Device B. (a) Bias spectroscopy at 
T = 6.5 K, i.e. above Tc ~ 1.7 K. Four-fold shell structure 
is observed with occupation numbers written in the large 
Coulomb diamonds corresponding to a full shell (zero occu- 
pation does not correspond to half- filling) . (b) Same gate 
region restricted to \Vsd\ < 3mV at T = 0.3 K showing in- 
elastic cotunneling features (blue arrows) as well as the onset 
of elastic quasiparticle cotunneling (horizontal green arrows). 
Every fourth diamond shows pronounced gate-dependent con- 
ductance peaks near \Vsd\ = 2A/e (vertical green arrow). 



possible Kondo ridge at this temperature in the normal 
state. The high Tc and critical field of the Nb-films pre- 
vented us from confirming the presence of a normal-state 
Kondo resonance, insofar as this resonance would already 
be suppressed by the magnetic field. The inelastic cotun- 
neling lines are seen more clearly in Fig. [3]^ a) which shows 
detailed measurements from the dashed rectangle in Fig. 
2(b). The cotunneling processes are depicted in Fig.jsjb- 
e), involving a weaker coupled orbital (orbital 1) [thin red 
line in Fig. |3jb)] and a stronger coupled orbital (orbital 
2) [thick blue line in Fig. |3jb)] split by 6. Such processes 
are allowed for all but the charge state corresponding to 
a filled shell [Fig. [sj^e)], consistent with Fig. [sj^a). As 
demonstrated by Holm et al.^^ a difference in tunnel- 
couplings to the two orbit als in the quantum dot, gives 
rise to a gate-dependence of the threshold for inelastic co- 
tunneling. With superconducting leads, such tunneling 
renormalization produces a gate-dependent shift of the 
unrenormalized threshold at eVgd ^ + 2A, as observed. 

In Fig. |4ja) we show conductance line-cuts through the 
center of each diamond 81-84 in Fig. [sj^a). The variation 
of the peak heights can be understood from the number 
of elastic (EL) and inelastic (INEL) cotunneling channels 
in each diamond, as summarized in Table [T[ There, the 
EL notation indicates whether the tunneling takes place 
through the higher- or lower-energy orbital. For example, 
in case of filling 1 [see Fig.jsj^b)], a total of 6 EL channels 
contribute: 2 from the upper and 4 from the lower orbital. 

In Fig. [4]^b) we show the results of a calculation of the 
lowest order nonlinear cotunneling conductance for the 
four different charge states in a single shell. The cal- 
culation involves the quasiparticle tunneling rates W^f 
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FIG. 3: (Color online) Device B. (a) Detailed bias spec- 
troscopy at T = 0.3 K of one four-electron shell [dashed blue 
rectangle in Fig. [2]^b)]. (b-e) Schematic energy diagrams il- 
lustrating inelastic cotunneling processes relevant for fillings 
of 1-4 electrons in a shell, e.g., corresponding to charge states 
81-84. The thick blue (thin red) level represents the strongly 
(weakly) coupled orbital, which both are coupled more weakly 
to the right electrode (drain) shown by a thicker barrier. 



FIG. 4: (Color online) Comparison of measured (a) and calcu- 
lated (b) conductance versus Vsd in the center of diamonds 81- 
84 in Fig.jsja). In the calculation 6 = 3A, ksT = 0.03 meV 
and addition energies [12.6, 12.6, 10.8 & 20.7] meV extracted 
from the width of the diamonds 81-84 have been used. 



between leads (3 = L^R and orbitals i,j = 1,2. For 
example, the rate for the process shown in Fig. jsjb) is 

p 4 

Wl^,^ = J^fpJ F^ (^)Ff + S)h {E) [1 -fniE^ S)] 

(1) 

with fa{E) = f{E - /ic.) the Fermi function, Tf{E) = 
rf\E - /i,|/v/(^-/il)2-A2, and Tf = nvp^A^ in 
terms of the tunneling amplitudes toc^i- Solving the 
steady-state rate equations we obtain the orbital occupa- 
tion numbers as a function of Vgd and the current is read- 
ily determined. In agreement with experiment one sees 
from Fig. |4] that the inelastic cotunneling peak is largest 
for diamond 82 due to the larger number of tunneling 
possibilities. Likewise the elastic 2A-peak is largest for 
diamond 83 due to the larger coupling of orbital 2 to the 
leads. The calculation does not reproduce the small shifts 
of some of the peaks seen in the experimental data, as 
well as the ratio between the amplitude of the elastic and 
inelastic peaks in diamond 81. We speculate that the lat- 
ter discrepancy will be removed by including higher order 
processes involving the strongly coupled orbital 2. The 
width of the peaks arise from a small smearing factor 
{t] = O.OlmeV) used in the DOS of the superconducting 
leads mimicking a constant inelastic scattering rate in the 



TABLE I: Number of cotunneling channels. 

Shell filling 12 3 4 

- il) ii) it) il) 

INEL 2 4 2 



contacts. We ascribe the larger width of the measured 
peaks to tunnel-broadening of the excited states which is 
not included in the calculation. 

For the results shown in Fig. Jl] we have used 
Ff = 0.07 meV, F^ 0.1 meV, Ff = 1.0 meV, and 
F2 = 1.8 meV, yielding an asymmetry factor in the order 
of rf / Ff 15 consistent with an upper bouncP^ 
of approximately 40 extracted from the Coulomb 
peak heights in the shell at 6.5K. The couplings are 
estimated from the sequential current at large positive 
(/+) and negative bias (/~) at the charge degeneracy 
point for adding the first electron in a four-fold de- 
generate shell. These currents are given by = 

rj-p-L / Rj~iL/ Ft /-j^R/ L .j^R/ L \ 

^ 2(r2 -^^rf/^+rf /^rf/^)+rf /^r2 11 n/i./ 0.0 n/i. 

in reasonable agreement with the experiment at the 
Coulomb resonance involving electron charge states 80 
and 81, itJp ^ 12nA/— 3nA. Moreover, the chosen 
couplings lead to a gate-voltage slope^ of the inelastic 
2A + ^ line of dS/d{ejVg) = 4 ^^(^2 - )/7r/7 ^ 0.084 
(average U = 12.0 meV), which agrees with the 
experimental result of 0.082 [average slope of the 
cotunneling lines in diamonds 81-83 in Fig. [sj^a)]. 
We note that since U A, the expression for the 
renormalization obtained in the normal state^ remains 
valid in the case of superconducting leads. Based on 
the couplings Tf used above, we estimate the Kondo 
energies (temperatures) in diamonds 81, and 83 to be 
^bTkm ^ 0.2 /ieV (2mK) and ksTx^ss ^ 0.03 meV 
(300 mK), i.e. much smaller than A ~ 0.28 meV, 
consistent with our observation of sub-gap structure 
rather than an enhanced zero-bias conductance pealP^. 

We now discuss measurements from a better coupled 
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FIG. 5: (Color online) Device B. (a) Bias spectroscopy at 
T = 0.3 K for a strongly coupled gate-region with pronounced 
change of the sub-gap structure in every fourth diamond. 
Numbers indicate the additional electron number in corre- 
spondence with Fig. [2] while the large/small diamonds are as- 
signed even (E)/odd (O) electron filling. The horizontal green 
arrows indicate the elastic quasiparticle cotunneling lines at 
eVsd = ±2A ^ ±0.5 meV together with the lower-lying lines 
corresponding to a single Andreev reflection, (b) Linear con- 
ductance (zero bias) versus gate voltage hinting that well and 
poorly coupled orbitals are filled consecutively, (c-d) Zoom at 
the gap structure around two Coulomb blockade resonances 
marked by arrows in (a) with the gate voltage in units of A. 



center of diamonds 3, 7, and 11. A similar negative 
differential conductance effect has been reported in a 
nanotube quantum dot with very different couplings to 
source, and drain electrodes^. 

In Figs. |5|c,d) we show a zoom-in of two neighbor- 
ing charge degeneracy points {N = 3/4 and N = 4/5) 
marked by vertical (blue) arrows in Fig.[5|a). Panel (d) 
shows a close resemblance to the behavior expect ed from 
MAR in the presence of a resonant level (F < /\> |19 | 24 | 25 | 
with the rounding of the sub-gap structure extending 
some 4A into the Coulomb blockade diamond. The data 
in panel (c), on the other hand, exhibit a qualitatively 
different sub-gap structure with pronounced negative dif- 
ferential conductance and a much stronger gate depen- 
dence which levels off at roughly 15-20A away from the 
charge-degeneracy point. 

We speculate that the unusual sub-gap features ob- 
served in odd occupied dots are caused by an inter- 
play between MAR and quantum Shiba states^^, present 
for spinful dots with ksTx < A. For strongly asym- 
metric couplings, these spin-induced bound states re- 
main pinned to the stronger coupled lead, at ener- 
gies inside the gap given roughly by the exchange 
coupling with this electrode. Therefore, new conduc- 
tance peaks away from the usual ±2A/n occur nat- 
urally in this scenario, and a bias-scan with the co- 
herence peaks of the weaker coupled lead gives rise 
to negative differential conductance at ±2A because 
of spectral weight transfer from the coherence peaks 
of the stronger coupled lead to the bound states. 



gate voltage region of device B, exhibiting a characteris- 
tic rounding of the elastic quasiparticle cotunneling and 
an unusual sub-gap conductance as seen in Fig.jsja). For 
charge states with three electrons in a shell, the sub-gap 
structure is especially pronounced and gate-dependent, 
indicating that this orbital is particularly well coupled to 
the leads. This is supported by the linear conductance 
data presented in Fig. [5|b), showing broad resonances 
where the effect is largest (diamonds 3,7 and 11). In 
diamond 3, the sub-gap conductance even exceeds 
the elastic cotunneling peak in the other diamonds 
at eVsd ~ with strong peaks at voltages 

which are different from the expected MAR positions 
at ±2A/n. Figure [sja) also displays marked negative 
differential conductance, seen as blue lines reflecting a 
local minimum in the current at eVsd ^ ±2A in the 



In summary, we have demonstrated how supercon- 
ducting electrodes lead to dramatic enhancement of 
cotunneling spectroscopy in carbon nanotube quan- 
tum dots. This revealed pronounced inelastic co- 
tunneling lines with marked gate-dependence caused 
by tunneling-induced level-shifts. Moreover, we dis- 
cussed the presence of negative differential conduc- 
tance and unusual sub-gap conductance in strongly cou- 
pled odd occupied diamonds. Further studies are re- 
quired to fully uncover the interesting interplay be- 
tween MAR and spin correlations in quantum dots. 
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